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Multifractal Wave Functions of a System with a Monofractal Energy Spectrum 

Masayuki Tashima* and Shuichi Tasaki 
Department of Applied Physics, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan 

We show the appearance of multifractal wave functions on a one-dimensional quasiperiodic system 
that has a monofractal energy spectrum. Using the Mantica technique, we construct the model as an 
inverse problem from the energy spectrum of a pure Cantor set. A relation between the critical state 
and the information dimension is proved and it is applied to the finite-size multifractal analysis. 
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1. INTRODUCTION 

Recently, the progress of experimental studies reveals interesting nature of the quasicrystal and 
new applications of the quasicrystal are expected to contribute to resolving a variety of problems. !) 
Under the circumstances of such experimental progress, theoretical study needs further development 
further to explain physical properties of the quasicrystal. Some previous studies have revealed inter- 
esting nature of quasicrystals and quasiperiodic systems as a theoretical model of the quasicrystal. 
Strange nature of wave functions is an especially attractive problem, for example a critical state and a 
multifractal wave function as an eigenstate. 2 ^ 

A fractal structure is also obtained in the energy spectrum of quasiperiodic systems. For example, 
the Fibonacci lattice is well known to have a Cantor-set-//&e spectrum (with zero Lebesgue measure 
and multifractality 3) ). It raises another question as to whether the energy spectrum is always fractal in 
quasiperiodic Hamiltonians and vice versa. This relation is an interesting problem, because the fractal 
dimensionality of the energy spectrum is connected to some physical properties of quasiperiodic sys- 
tems, for instance, the temporal autocorrelation function. 4 ^ Some study showed the above relation in a 
particular model, 5) but no general proof is available. 

One of the difficulties of theoretical study may be caused by the structure of quasiperiodic sys- 
tems. A quasiperiodic function is designed to be a function that can be uniformly approximated by a 
Fourier series with a finite number of pairwise incommensurate base frequencies. 6 ^ Because of the in- 
commensurability, quasiperiodic functions never have the translational periodicity and need a special 
construction technique. Well known ways include the inflation rule, 3) an incommensurate potential, 
and the projection method. 7 ^ They were used in many previous works, for instance the Fibonacci lat- 
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tice, the Thue-Morse lattice, the Harper model (ID), and the Penrose tiling (2D). They succeeded 
in explaining some aspects of the quasiperiodic system. However, the absence of simple symmetries 
of crystals makes further study of quasiperiodic systems difficult, and a new approach to studying 
quasiperiodic system is strongly needed. 

Mantica proposed a completely different way, an inverse-problem approach. 8)9) The Mantica tech- 
nique is a way of constructing a quasiperiodic system from a multifractal spectrum. It offers us another 
way of studying quasiperiodic systems. Gurneri and Mantica used the technique and revealed the re- 
lation between the fractal energy spectrum and the anomalous diffusion of wave packets. 10) 

In the present paper, we construct a one-dimensional quasiperiodic system with a raonofractal 
energy spectrum by the Mantica technique and study the mwWfractality of its wave functions. The 
Mantica technique enables us to obtain an eigenstate of the singular continuous spectrum more easily 
than ordinary ways because we first set the energy spectrum and then calculate the corresponding 
Hamiltonian. We choose a pure Cantor set (with zero Lebesgue measure and monofractality) as the 
energy spectrum. We conjecture that the dimension of the monofractal spectrum is related to the fractal 
dimension of the wave functions. Moreover, we prove the relation between the information dimension 
and the index of the critical wave function. Using the relation, we evaluate the multifractal spectrum 
for a finite number of lattice systems. 

2. MANTICA TECHNIQUE 

We construct a semi-infinite one-dimensional quasiperiodic system as an inverse problem through 
the Mantica technique. It enables us to construct quasiperiodic systems from some kinds of multifrac- 
tal spectra. Before presenting our results, we briefly review the Mantica technique introduced in the 
1990s. 8)9 ) This is a technique of constructing a tridiagonal Hamiltonian possessing a given singular 
continuous spectrum belonging to the class of iterated function systems (IFS). 

The Schrodinger equation for a semi-infinite one-dimensional tight-binding model is generally 
given by 



where if/„(E) is the wave function at site n with the energy E, e„ is the on-site potential, and t„ is the 
hopping element between the sites n - 1 and n. This equation constitutes a recurrence formula, so that 
if/ n for energy n can be calculated from deciding the first two terms and all coefficients [e n , t n \. 



On the one hand, a set of orthogonal polynomials {P n (x)}™ with respect to a measure n(x) (//(O) = 
and Ml) = 1), 



E\p n (E) = t n+ n// n+ i(E) + e n i/f n (E) + t n i// n -i(E), 



(1) 




(2) 
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satisfies the same form of equation 1 1) 

xP n (x) - lU) + a n P n (x) + b n P n -i(x), (3) 

where 

a n = \ xP n (x) 2 dfi(x) 
Jo 

and 



b n = | xP n -i(x)P n (x)dn(x). 
Jo 



H) 

We can find the following correspondence by comparing Eq. (1) to Eq. (3): 

"An <-> 

E <-» X 

H(E) <-» fi(x) 
e n <-» a n 
t n <-> b n , 

where //(fi) is the integrated density of states. Consequently, solving Eq. (1) is replaced by solving the 
three-term recurrence formula (3). 

In general, computing a whole set of orthogonal polynomials {P n } associated with the measure 
is not an easy endeavor, but with the Mantica technique we can calculate the polynomials {P n } by 
using the invariability of IFS under an appropriate affine transformation 

x -> 5x + /?, (4) 

where 5 and fi are real constants. The integral with a measure fi(x) satisfies 

f(x)dfi(x) = ^ TTi JT f(5ix+pi)dn(x), (5) 

where the weights Jij > with £,-7T; = 1 for any continuous function /. The parameters 6j, Pi, and ni 
are determined from a preassigned spectral measure fi(x) that belongs to IFS. In the Mantica technique, 
the parameters a n (- e n ) and b n (= t n ) can be calculated through Eq. (5) and the polynomials P n are 
also calculated from the initialization P_i = and Pq = 1. 8) We can thereby obtain the wave function 
if/ n (E) for an integrated density of states fj.(E) that belongs to the class of IFS. 

We calculate the first-order polynomial P\(x) as an example. The first-order polynomial P\(x) 
must has the form 

P l (x) = C l x + C , (6) 
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where Co and C\ are real constants. The orthogonal relations for the first polynomial are given by 

f Po{x)P x {x)dn(x) = Ci{x) + CQ = (7) 
Jo 

and 

f Pi(x)Pi(x)dn(x) = Ci 2 (x 2 ) + 2CiCo(x> + Co 2 - 1. (8) 
Jo 

Here (x n ) is 

<x"> = f x n d f i(x). (9) 



Thus the constants are given by 



and 



Ci - 1 (10) 

V<* 2 > - (x) 2 



Co = C!<X>. (11) 

We can calculate the constants by determining (x 2 > and (x). They are obtained by the affine transfor- 
mation Eq. (5) as follows: 

<x>= f xdn(x) = Ym[ (Six+ftdnix)) = Ymidiixy+fc). (12) 
Jo Y Jo i 

Therefore (x) is described by the parameters tt,-,/?,, and 8i 

< x > = i v P ^ 13 > 
1 -ZiXiSi 



We can determine (x 2 > in the same way: 

As a result, the first-order polynomial Pi(x) = C\x + Co is obtained from the parameters /r ( , /?,-, and 5,-. 
The higher-order polynomials P n (x) (= if/ n (E)) are also computed in the same manner. 

In our study we choose as the energy spectrum the Devil's staircase (Fig. 1), which is described 
as IFS in the form, 



yu(x) = 



\nOx) (0 < x < i), 



3 
2 

2 v 3 - ■* ^ 3 



a < x < h, d5) 



\nOx - 2) + 3 (§ < x < 1), 
where x is the energy normalized in [0,1]. This corresponds to a pure Cantor set (middle third re- 
moved). This is a monofractal, since the generalized dimension D q of the measure [i is D q - Dq - 
log 2/ log 3 0.631. The parameters for the affine transformation (5) are derived from the spectrum 
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Fig. 1. The Devil's staircase. 



(15) as 

r l r l/3 f 2/3 r l 

/(x)JyU(x) - /(X)^(X) + f(xW( X ) f(x)dn(x) 

JO Jo Jl/3 J2/3 

■jC(5'(I)*H^))** 

Thus, the parameters are 6\ = 82 = 1/3, n\ - tt2 = 1 /2, jSi - 0, and = 2/3. 

3. QUASIPERIODICITY OF THE TIGHT-BINDING MODEL 

In this section, we construct a one-dimensional system from the Devil's staircase and confirm its 
quasiperiodicity. The semi-infinite lattice system can be obtained from the recursive sequence of the 
Mantica technique, but calculating a very large number of polynomials corresponding to the number 
of lattice sites is a hard work for computer. Hence we choose 2 13 = 8192 as the site number. Setting 
the site number finite never changes the boundary condition. ll) 

First we investigate the power spectrum of the on-site potential e n and the hopping element t n in 
order to examine the quasiperiodicity of the system. First 50 data points are excluded so as to eliminate 
the boundary effect at n - 0. The Fourier spectrum of e n and t n are given by 

|F e | 2 = |£(e„ - (e))e ink \ 2 (17) 



and 



\F t \ 2 = \Yj,t n -(t))e ink 



(18) 



where (e) and (t) are the arithmetic means of e n and t n , respectively (see Fig. 3 and Fig. 3). The 
spiky structure indicates a long range order. The positions of the peaks without a simple proportional 
relation, on the other hand, indicates that it has no periodic arrangement. We thereby claim that the 
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Fig. 2. The power spectra of the (a) on-site potential e„ and (b) the hopping element t n . The power spectrum 
of hopping element t„. Since peaks are at k/n = 0.5741 . . . , 0.7745 . . . , 0.3404 . . . , the Fourier spectrum is 
not described by a simple proportional relation. A clear hierarchical structure is also observed. 



system is quasiperiodic. ' 

The potential and the hopping element of our system do not take only two values in real space, 
while the Fibonacci lattice and the Thue-Morse lattice consist of two elements as the potential or the 
hopping (the quasiperiodicity only appears in its arrangement). The similarity between our system 
and these quasiperiodic systems is unfortunately not observed from the viewpoint of Hamiltonian 
elements. Moreover, the potential of our model is not described by a cosine curve with an irrational 
number, such as in the Harper model. 

4. MULTIFRACTAL ANALYSIS OF CRITICAL STATE 

In present section, we describe how we can characterize the fractality of the wave functions of 
our system. Wave functions are generally classified into three types, namely extended, localized, and 
critical. 14) An extended state is defined as 

f W)U?~L D , (19) 

J\r\<L 

where L is the system size and D is the spatial dimension. It is similar to the Bloch state. A localized 
state is 

f Wftfdf-I?. (20) 

J\r\<L 

The last type of the wave function, the critical state is neither localized nor extended: 

r w?fdf~L\ (2i) 

J\r\<L 

where the index v is in the range < v < D. A typical wave function in the critical state may be a 
power-law type function if/(r) ~ \r\ y with v < D/2 or a self-similar function. 14) 
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Equation (21) can be rewritten for our one-dimensional discrete system as 

N 

s(N) = Yj^ 2 ~ NVn > < 22 ) 

where N is the site number and is a value of the index v for the A-site system. The logarithm of Eq. 
(22) leads 

In 5 (AO ~ v N \nN. (23) 

We can estimate the index in the range (0, 1) as a gradient in a plot of lnS(A0 against In N. 

Another way of characterizing the critical state is the multifractal analysis. The multifractality 
is represented by the singularity index, or the Lipschitz-Holder exponent, a q , and the multifractal 
spectrum f{a q ). l ^ l6) The singularity index a q describes the local degree of singularity (local fractal 
dimension) and the multifractal spectrum f(a q ) is the fractal dimension of the support which has the 
singularity a q . In other words, the curve f(a q ) against a q means a relation between the local fractal 
dimension characterizing the system and the global fractal dimension of the spatial distribution of the 
singularity. 

The multifractal analysis has been applied to study of the nature of wave functions. The analysis 
shows that an extended state has f(a q - 1) = 1, a localized state has f(a q = 0) = and f(a q - oo) = 
1, and a critical state only has a smooth convex multifractal spectrum f(a q ) in the range [c^min, a max ]. 
For example, the multifractal spectrum for a Fibonacci lattice was obtained by T. Fujiwara et al. 2) They 
obtained an asymmetric shape of f(a q ) and the limits f(a max = ar_oo) - and /(cK m i n - a +co ) + for 
the wave functions at the edge of the energy spectrum. The multifractal spectra f(a q ) of the critical 
states of other quasiperiodic systems are also smooth. 

In the present paper, we consider a discrete lattice system, and therefore the multifractal formalism 
is modified to the following equations. The spatial distribution <2, (which is described as the probability 
measure in mathematics) at cell i with size / is given by 

/n lS M \p k (x)\ 2 

& 77 = tl . , (24) 



where the summation in the denominator is for normalization and N is the number of the sites on the 
whole lattice. The participation function Z q is defined by 



7 \ m 

H 



(} 4 



(25) 



where m = N/l is the number of the cells. In Eq. (25), different parts of the distribution Qi can be 
stressed by the parameter q. If q has a large positive value, the participation function Z q is dominated 
by the regions corresponding to larger values of <2,-, whereas if q is a negative large number, Z q is 



7/14 



J. Phys. Soc. Jpn. Full Paper 



D q = - — 7 hm lim — ^— — . (26) 



dominated by contributions from small values of Q,-. The generalized dimension D ? is given by 

1 lnZ ? (//AT> 
lim lim 

q - 1 /v^c» <r->e mill In e 

A lattice system has the lattice constant as the minimal value, and therefore the extremum e has to be 
fixed at the lower cut off length scale e m ; n = l/N = 2 IN for our system. The singularity index a q and 
the multifractal spectrum f(a q ) are calculated from the D q as 

a « = ^[ ( *- 1)D J (27) 

and 

f{a q ) = qa q - (q - \)D q . (28) 

As is easily seen from Eq. (27) and Eq. (28), the multifractal spectrum for a monofractal reduces to 
the dimension of the monofractal Df = f(a q ) = a q . 

If q = 1, the generalized dimension D q is called the information dimension D\ 

i ^ tiv // 



Di = lim > - In Of — . ( 29 ) 

The information dimension Di coincides with the index of the critical state given in (22) (see 
Appendix). The information dimension for the finite number of lattice system Df is described as 

Df = v N + 6D { (N), where lim SD^N) = 0. (30) 

N— >co 

We can evaluate the multifractal spectrum for the finite number of lattice system through the n depen- 
dence of Df - vn. 

5. RESULTS AND DISCUSSION 

A typical wave function is illustrated in Fig. 3. It clearly lacks periodicity. To characterize the 
structure of the wave function, we study the critical state and the multifractality. First 50 data points 
are also excluded from both analysis. 

Let us first describe the analysis given in Eq. (23). We evaluate S(N) for the wave function illus- 
trated in Fig. 3. The resulting plot in Fig. 4 yields the index vn - 0.886. It indicates that this wave 
function is in the critical state. 

Wave functions in the critical state may be a power-law type function or a self-similar function. ' 
The wave function illustrated in Fig. 3 might be a series of power-law-type localized functions. In 
order to check the possibility, we divided the wave function into five clusters, where a cluster is 
defined as a region of large amplitudes between sites of small amplitudes as indicated in Fig. 3 and 
applied the power-law approximation to each cluster. However, the clusters were not well described 
by the approximation. In addition, the function in Fig. 3 is not a simple self-similar function, that 
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Fig. 3. Eigenfunction for the energy x = 1/3. There are five clusters except for the clusters around the bound- 
ary sites n = and n = 2 n = 8192. 




In[nl 

Fig. 4. S (/V) for the eigenfunction for the energy x = 1/3. 

was reported on a Fibonacci latticed Hence, we analyzed them from the other viewpoint, namely the 
multifractal analysis. 

We then investigate the multifractality of the wave functions. Figure 5 shows the multifractal 
spectrum of the wave functions of the system. We thereby conclude that the wave function of the 
system with a monofractal energy spectrum exhibits multifractality. The maximum value of f(a q ), the 
dimension of the support, is unity. The equivalence of the spatial dimension and the dimension of the 
support means that there is no areas with zero amplitude. We calculate the average of the limitation 
values of the multifractal spectrum as /(a m in) - 0.038 at a m i n 0.58 and f(a max ) - at a max 2.7. 
The multifractal spectrum f(a q ) exhibits the maximum value unity at a q - \ .2. Any simple relations 
between the dimension of the energy spectrum D cantor = log 2/ log 3 0.631 and the multifractal 
spectrum f(a q ) are not observed in Fig. 5. 

The multifractal analysis for the finite-size lattice system is needed to be treated with discretion. 2) 
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Fig. 5. The multifractal spectra f(a q ) are plotted for the wave functions with the energies x = 1/27 (solid 
line), x = 1/9 (broken line), x = 2/9 (long dashed line), x = 1/3 (dotted line), and x = 2/3 (dash-dotted 
line). 
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Fig. 6. The declination Dj - v N against the number of site n for the energy x = 1/3. 



Thus we confirmed the convergence of the multifractal spectrum f(a q ) by Eq. (30). We can see the 
convergence of D 1 ^ - vn against /V in Fig. 6. 

6. SUMMARY 

We constructed a quasiperiodic system from the pure Cantor set spectrum by the Mantica tech- 
nique and proved the relation between the critical state and the information dimension. The preas- 
signed spectrum is monofractal. Nevertheless, the wave functions of our system exhibits the multi- 
fractality. In this work, we have not observed a simple relation between the dimension of the energy 
spectrum and the multifractal spectrum of the wave functions. More study is needed to understand the 
fractal nature of wave functions of quasiperiodic systems. 
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Appendix: A RELATION BETWEEN THE CRITICAL STATE AND THE INFORMATION 
DIMENSION 

The critical state for a finite number of lattice sites N is defined as 

N 

J]m 2 = C N N VN , (A-l) 
(=1 

where i/>; is the wave function at site i, Cn is a positive constant and, vjv is the index for Af-sites lattice 
system. We suppose that vn coincides with v in the limiting case 

lim vn - v. (A-2) 

N->co 

We divide the whole system into m parts with size / - N/m and assume that the summation of the 
wave functions of z'th cell Ni satisfies 

Ni 

^W 2 = Cif», (A-3) 
j 

where C; is a positive constant and v,- is a critical index of the z'th cell. The summation of Eq. (A-3) 
over m parts should reduce to Eq. (A-l): 

m Nj m 

Z Z ^ |2 = Z Q r ' = CnNVN - (A4) 

i 7=1 i=l 

The spatial distribution Q { at the ith cell is given by 

G, Urz*itf-cvjv» (A5) 



and Qi is normalized as 



y Qi [L l — I — Vc < r» = i. (a-6) 

j-( \Nj C N N V » j-( 
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(A-7) 



The summation in the above equation is rewritten by 



1 m 1 m 

—— V d l v - In (Cj f) - — — In (Cjy A™) V Cj Z v ' . 



(A-8) 



i=l 



The second term of Eq. (A-8) reduces to 

— ln(CiviV v,v )2]c«Z y ' ^ln^iV^) -lnC^ + v^vlniV. 



i=i 



Therefore the information dimension Di is given by 
1 



(A-9) 



Di = lim 



w->oo In N - In Z 

In C^ + y^ In Z 

= y + lim 



In + In N - g Cj l v ' In C« + In Z g v/C,- Z v 



1 



l£i ' 



(A- 10) 



In d + In Z ^ Vid l v 

(=1 



#-»«> lnA^-lnZ lnA^-lnZC WJ 
We now define K] max as the maximum value of {X,} in order to evaluate a behavior of the absolute 
value |lnC,-|: 



max{|lnCj|} = [|lnCj|] 

max- 

i 

We then obtain the following inequality: 

m m 

£ Cj Z v ' In Cj < ^ Q Z v < [|ln C«|] max - [|ln Cj|] max C w ZV V ~ 



i=i 



(A- 11) 



(A- 12) 



Then then the third term of Eq. (A- 10) converges to zero: 



1 



1 



lim 

n^oo In N- In IC N N VN 



£cj/ v -lnCj 



i=i 



< lim 



1 



1 



a^oo In N - In ZC^v A"* 

[|lnCj|] max 
= hm ; — rrr: — ; — — - 0. 



[|lnCj|] raax OvA^ 



(A-13) 



In [N] - In [Z] 

The forth term in Eq. (A- 10) is evaluated in the same manner. It also converges to zero: 

lim — ^ L_ f ViCl r- < hm = o. (A-14) 

a^co In Ar- In IC N N VN j-{ JV->°° In A^ - In Z 

We thus arrive at the relation between the information dimension D\ and the index v: 

Di(=ai =/(a0) = v. (A-15) 

If a state is localized, Di = v = and if a state is extended, D\ - v - 1. These results are consistent 
with a previous research. 14) We can obtain the same relation for a continuous system in the same way. 
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= v N + 6Di(N) 



1 1 ( ^ N VN ) 

VN + rrr-n t^~^ \\c N \nc N cr in q 

mN - In/ CnN Vn yj-f n 

) v N C N ViCil ' . 
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(A-16) 



lnN-lnlC N N v 

The each term in the brackets can be interpreted as representing a roughness of a local structure. 
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